Constraints on Parton Charge Symmetry and Implications for Neutrino Reactions 
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For the first time, charge symmetry breaking terms in parton distribution functions have been 
inciuded in a giobai fit to high energy data. We review the results obtained for both valence and sea 
quark charge symmetry violation, and we compare these results with the most stringent experimental 
upper limits on charge symmetry violation for parton distribution functions, and with theoretical 
estimates of charge symmetry violation. The limits allowed in the global fit would tolerate a rather 
large violation of charge symmetry. We discuss the implications of this for the extraction of the 
Weinberg angle in neutrino DIS by the NuTeV collaboration. 
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I. INTRODUCTION 

Charge symmetry is a restricted form of isospin invari- 
ance involving a rotation of 180° about the "2" axis in 
isospin space. For parton distributions, charge symmetry 
involves interchanging up and down quarks while simul- 
taneously interchanging protons and neutrons. In nuclear 
physics, charge symmetry is generally obeyed at the level 
of a fraction of a percent Q, Q • Charge symmetry viola- 
tion in parton distribution functions (PDFs) arises from 
two sources; from the difference 5m = — m u between 
down and up current quark masses, and from electromag- 
netic (EM) effects. 

Since charge symmetry is so well satisfied at lower en- 
ergies, it is natural to assume that it holds for parton 
distributions. Introducing charge symmetry reduces by 
a factor of two the number of PDFs necessary to describe 
high-energy data. In addition, as we shall see, there is 
no direct experimental evidence that points to a substan- 
tial violation of charge symmetry in parton distributions. 
For this reason, until recently all phenomcnological PDFs 
have assumed charge symmetry at the outset. In Sec. 
ITT1 we review the experimental evidence for charge sym- 
metry in PDFs. Recent experiments allow us to place 
reasonably strong upper limits on parton charge sym- 
metry violation. In Sec. IIIII we review some theoretical 
estimates for parton charge symmetry violation (CSV) 
for valence quarks, and we compare this with the exper- 
imental limits. In Sec. IIVI we review the recent global 
fit by Martin et ah, Q that includes for the first time 
the possibility of charge symmetry violating PDFs both 
for valence and for sea quarks. We compare the results 



of this phenomcnological fit with the experimental limits 
on parton CSV as well as with theoretical estimates of 
both the magnitude and sign of charge symmetry violat- 
ing parton distributions. 

In Sect. we review the effect of isospin violating 
PDFs on neutrino DIS measurements. In particular, we 
concentrate on the extraction of the Weinberg angle by 
the NuTeV collaboration p|. We show the magnitude 
of the effects predicted by theoretical calculations, and 
we contrast this with the magnitude of effects allowed 
by the phenomenological CSV PDFs extracted by the 
MRST group 0. 



II. EXPERIMENTAL LIMITS ON PARTON 
CHARGE SYMMETRY VIOLATION 

There are no direct measurements that reveal the pres- 
ence of charge symmetry violation in parton distribution 
functions. At present, we have only upper limits on the 
magnitude of charge symmetry violation. 

The most stringent test of parton charge symmetry 
comes from comparing the structure function F% mea- 
sured in neutrino induced charged current reactions, and 
the structure function for charged lepton DIS, both 
measured on isoscalar targets N - In leading order, as- 
suming parton charge symmetry, the structure functions 
have the form [§J 



F^ ±N °(x) 
Q{x) 



^> + | (c(x) + c(x) - s(x) - S(x)) 

Q{x) ± x (s(x) — s(x) + c(x) — c(x)) 

x(q j (x)+q j (x)) (1) 

j—u,d,s,c 
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In Eq. ^ FY N ° 1S the structure function for charged- 
current processes induced by neutrinos (antineutrinos) 
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on an isoscalar target. Charge symmetry violation and 
NLO effects are not included in this equation. In the 
limit of exact charge symmetry, there is a simple rela- 
tion between the charged-lepton and neutrino structure 
functions, corrected for heavy quark effects. This rela- 
tion is defined as the "charge ratio" R c (x,Q 2 ) or, as it 
is sometimes termed, the "5/18" 1 rule." The factor 5/18 
is simply understood as the average of the squares of the 
light quark charges, relative to the weak charges. Ex- 
panding R c to lowest order in the (presumably small) 
charge symmetry violating terms gives 



Rc(x) 



F2 No (x) + x (s(x) + s(x) - c{x) - c(x)) /6 



m WNo (x)/l8 



3x (Su(x) + 5u(x) — dd(x) — 5d(x)j 
1+ 10Q(x) 



Eq. [5] introduces the CSV parton distributions, 

Su(x) = u p (x)-d n {x); 
Sd(x) = d p {x)-u n {x), 



(3) 



with analogous relations for antiquarks. The quantity 
R c {x) in Eq. is defined using F% °, the average of 
neutrino and antineutrino, F 2 , charge-changing struc- 
ture functions; it also requires knowledge of strange and 
charm PDFs. Deviation of R c (x) from unity would be 
an indication of a non-zero charge symmetry violating 
contribution. 

The most precise neutrino measurements were ob- 
tained by the CCFR group ,6j, who extracted the F 2 
structure function for neutrino and antineutrino inter- 
actions on iron using the Quadrupole Triplet Beam at 
FNAL. This can be compared with several measurements 
of the F 2 structure functions from DIS reactions using 
high-energy muons or electrons. The most precise mea- 
surements were obtained by the NMC group 0, El , who 
measured F 2 structure functions for muon interactions on 
deuterium at muon energies = 90 and 280 GeV. Ear- 
lier measurements were obtained by the BCDMS muon 
scattering experiments on deuterium [t| and carbon [lOj , 
and electron scattering results from SLAC [H[l3. 

Precision measurements of the charge ratio require a 
significant number of corrections. It is necessary to know 
the relative normalization between charged-lepton and 
neutrino cross sections. Heavy quark threshold effects 
are very important, as the most important neutrino en- 
ergies are sufficiently small that one must correct for the 
finite mass of the charmed quark EH- The CCFR neu- 
trino cross sections were taken on iron targets, so one 
must correct for nuclear effects (Fermi motion at large x, 
EMC effects at intermediate x and shadowing at small 
x El E E OH El). There are additional correc- 
tions since iron is not an isoscalar target. Eq. also 
requires contributions from strange and charmed quarks. 
The CCFR [U and NuTeV |2l| groups have extracted 
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FIG. 1: The charge ratio of Eq. Hvs. x using CCFR neu- 
trino data of Ref. [g combined with muon DIS data. Solid 
circles: NMC fj, — D data, Ref. Q; open triangles: muon 
measurements from the BCDMS group 011 D, Ref. i and C, 
; solid triangles: SLAC electron scattering data Refs. 



Ref. JT3; 



strange (antistrange) quark distributions from the cross 
sections for opposite sign dimuons in reactions induced 
by neutrinos (antineutrinos). 

Fig-Dplots the charge ratio R c of Eq.|2|vs. x. The cir- 
cles are the NMC/CCFR ratio. The open triangles are 
the BCDMS/CCFR charge ratio, where BCDMS repre- 
sents the muon scattering results of the BCDMS group on 
deuterium |9j and carbon [10j . The solid triangles are the 
SLAC/CCFR charge ratio, where SLAC denotes electron 
scattering results of the SLAC group 0, . Since the 
CCFR measurements were obtained from neutrino-iron 
scattering, and the NMC measurements were from fi — D 
reactions, the NMC measurements were "converted" to 
the equivalent reactions on iron by multipl ying by the 
ratio of fi — Fe to fj, — D cross sections [23, |23| . In ad- 
dition, a correction for strange quarks and charm quarks 
was applied to the numerator of Eq. [21 

In the region 0.1 < x < 0.4 the charge ratio test is 
consistent with unity, with errors in the range 2 — 3%. 
Solving Eq. |21 for the CSV parton distributions, 



[Su(x) + Su(x) — 5d(x) — 5d(x)) 



(4) 



one can set an upper limit to parton CSV effects in this x 
range at about the 6 — 9% level. For larger values of x the 
upper limit on CSV effects is substantially greater. This 
is due both to poorer statistics and to the large Fermi 
motion corrections needed for the heavy target at large 
x. 

At values x < 0.1, the charge ratio appeared to devi- 
ate significantly from unity. The discrepancy increased 
monotonically as x decreased, and approached 15% at 
the smallest values of x j24|- Several suggestions were 
advanced to explain this discrepancy, including anoma- 
lously large contributions from strange quarks [25|, and 
surprisingly large contributions from sea quark charge 
symmetry violation 26]. Eventually the CCFR group re- 
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analyzed the neutrino data, and after the new analysis 
the small- a: discrepancy disappeared |2jj , largely as a re- 
sult of two factors. The first was the treatment of charm 
mass corrections. In the initial analysis, these were in- 
corporated using the "slow rescaling" hypothesis due to 
Georgi and Politzer j^. The re-analysis involved NLO 
calculations, which differed significantly from the slow- 
rescaling procedure at small x - as had been suggested 
by Boros et al. [2j|. Since the NLO effects depend on 
both x and Q 2 , the final results cannot be plotted simply 
versus Bjorken x. 

The second significant effect involved the separation of 
the F 2 and F 3 structure functions in charged-current neu- 
trino DIS. The sum of neutrino and antineutrino charged- 
current DIS cross sections contains a linear combination 
of neutrino F 2 and F 3 structure functions, 



2„v 



d 2 <7 



1„v 



d 2 o 



dxdy dxdy 



~ 2(1 -y-y 2 / '2) F 2 (x,Q 2 ) 

+ (y-y 2 /2)AxF 3 (x,Q 2 ) . (5) 



In Eq. F 2 is the average of the F 2 structure functions 
for neutrinos and antineutrinos, and in leading order (as- 
suming charge symmetry), AxF 3 = 2x(s + s — c — c). The 
structure functions are multiplied by coefficients that de- 
pend on the invariant y = p-q/p-k, where k(p) is the four 
momentum of the initial lepton (nucleon), and q is the 
four momentum of the virtual W exchanged in the inter- 
action. In Eq. [SJwe have dropped terms of order m 2 N /s, 
and for simplicity the longitudinal/transverse ratio R has 
been set to zero. 

In the initial data analysis, the data for a given x bin 
was averaged over all y, and the AxF 3 structure func- 
tion was estimated using phenomenological PDFs. The 
re- analysis included the x, y and Q 2 dependence of the 
cross sections. In this way the group was able to ex- 
tract both F 2 and AxF 3 . The experimental values for 
AxF 3 differed substantially from the phenomenological 
predictions. This affected the extracted values for the F 2 
neutrino structure functions. The description of charm 
production also plays a si gnifi cant role in determining 
the value of AxF 3 (HSIIJli. The combined effect of 
the NLO treatment of charm production, and the model- 
independent extraction of AxF 3 removed the small-x dis- 
crepancy. The charge ratio R c of Eq. is now unity to 
within experimental errors, even at small x. Since the 
NLO treatment depends separately on Q 2 and x, it is no 
longer possible to display the results in a Q 2 -independcnt 
plot such as Fig. ^ 



III. THEORETICAL ESTIMATES OF PARTON 
CHARGE SYMMETRY VIOLATION 

Both the quark mass difference and EM effects re- 
sponsible for parton charge symmetry violation represent 
small changes to quark PDFs. Several theoretical esti- 
mates of parton CSV for valence quarks can be obtained 



by examining simple models for parton distribution func- 
tions, then observing how these change upon applying the 
operations of charge symmetry. 

The Adelaide group H, |33, Hj, |35j developed a method 
for calculating twist-two valence parton distributions 
through the relation 

<? v (:rV) = M \(X\i; + (0)\N)\ 2 6(M(l - x) - p+) . 
x 

(6) 

In Eq. -0+ = (1 + 03)^/2 is the operator that re- 
moves a quark or adds an antiquark to the nucleon state 
\N), [x 2 is the starting scale (where QCD is best ap- 
proximated by a valence dominated quark model) for 
the quark distribution, \X) represents all possible fi- 
nal states that can be reached with this operator (i.e., 
\X) = 2q, 3q + q, Aq + 2q, . . .), and p\ is the plus compo- 
nent (p + = p 3 + E(p)) of the momentum of the residual 
system. We note that the advantage of using Eq.0is that 
the correct support for the PDF is assured, regardless of 
the approximation used to evaluate the nonperturbative 
matrix element (X\ip+(0)\N). 

Sather j3(J assumed that the dominant result of elec- 
tromagnetic effects was manifested in the neutron-proton 
mass difference SM = M n — M p . At the low Q 2 appropri- 
ate to quark models, the largest contribution to valence 
PDFs at large x arises from configurations where one 
quark is removed from three valence quarks, leaving a 
residual diquark. Since CSV corresponds to interchang- 
ing up and down quark, and neutron-proton labels, va- 
lence quark CSV is obtained by examining how Eq. [B] 
changes with 8m and SM. 

Restricting Eq. to the residual diquark state X = 2 
(valid for large x) and neglecting the dependence on 
transverse momentum in the 5-function, Sather obtained 
analytic relations between valence quark CSV and deriva- 
tives of the valence PDFs, 



Sd v (x) = g?^ (x) — it" (x) 

= -™±\xd (x)] 
~ M dx [ v[ )l 

Su v (x) = Uy(x)—d%(x) 
SM ( d 
~~M 



5m d 
M dx 



(ly(x) 



, L xu v (a;)] + — Uy(x) 
dx dx 



(J) 



In Eq. 0, M is the average nucleon mass, SM = 1.3 
MeV is the n-p mass difference, and Sm = m<i — m u ~ 4 
MeV is the down- up quark mass difference. From Eq. 
\7\ at large x the down quark valence distribution in the 
proton is expected to be larger than the up quark distri- 
bution in the neutron; similarly, the down quark distri- 
bution in the neutron is larger than the up quark distri- 
bution in the proton at large x. The lowest moment of 
the valence quark distributions is fixed by quark normal- 
ization, since 



(8) 



dxSd v (x) = / dxSu v (x) = . 
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FIG. 2: Valence quark CSV contributions, xSq v (x) vs. x. 
Solid line: xSu v ; dash-dot line: xSd v . Calculated by Rodi- 
onov et al, Ref. |22l using MIT bag model wavefunctions, 
evolved to Q 2 = 10 GeV 2 . 

A violation of Eq. [S] would be equivalent to changing 
the total number of valence quarks in the proton and/or 
neutron. Eq. 00 shows that if 5d v (x) is positive at large 
x, it must therefore be negative at small x such that the 
integral over all x vanishes; similarly, the quantity Su v (x) 
must also change sign at small x. 

Sather's results do not include effects from configura- 
tions with more than two quarks in the final state. In ad- 
dition, his analytic results include only the contribution 
from quark longitudinal momentum (effects from quark 
transverse momentum in Eq.[S]are neglected). Rodionov 
et al., 37] explicitly included the effect quark transverse 
momentum within the MIT bag model. They did not use 
the approximate equations of Sather, but included quark 
and nucleon mass differences directly in Eq. |SJ Their re- 
sults for valence quark CSV PDFs are shown in Fig. [2] 
The solid curve is x5u v (x), while the dot-dashed curve 
is xSd v (x), both evolved to Q 2 = 10 GeV 2 . Qualita- 
tively, the results of Rodionov et al. are very similar to 
Sather's. The sign and relative magnitude of both Sd v (x) 
and Su v (x) are quite similar in both calculations, and the 
second moment of both distributions is equal to better 
than 10%. 



IV. PHENOMENOLOGICAL CHARGE 
SYMMETRY VIOLATING PDFS 

Because CSV effects are typically very small at nu- 
clear physics energy scales 0, , from the lack of direct 
evidence for violation of parton charge symmetry, and 
because theoretical estimates put parton CSV at below 
the 1% level @, all previous phenomenological parton 
distribution functions have assumed the validity of par- 



ton charge symmetry at the outset. However, Martin, 
Roberts, Stirling and Thorne (MRST) have recently 
studied the uncertainties in parton distributions arising 
from a number of factors, including isospin violation. 

The MRST group chose a specific model for valence 
quark charge symmetry violating PDFs. They con- 
structed a function that automatically satisfied the quark 
normalization condition of Eq. [SI namely: 

5u v (x) = —5d v (x) = nf(x) 
f{x) = (1 - x^x- - 5 (x - .0909) . (9) 

The function f(x) was chosen so that at both small and 
large x, f(x) has the same form as the MRST valence 
quark distributions, and the first moment of f(x) is zero. 
The functional form of the valence CSV distributions 
guaranteed that Su v and Sd v would have opposite signs at 
large x, in agreement with the theoretical results shown 
in Fig. H 

Inclusion of a valence quark CSV term changes the 
up and down quark distributions in the neutron from 
their charge symmetric counterparts in the proton. This 
changes the momentum carried by valence quarks in the 
neutron from those in the proton, since the total mo- 
mentum carried by valence quarks is given by the second 
moment of the distribution, e.g., the momentum carried 
by up valence quarks in the neutron 

Uy = / XU"(X) dX . 

Jo 

Because the total momentum carried by valence (up plus 
down) quarks in the neutron, U™ + D™, is rather pre- 
cisely determined, MRST chose a functional form that 
kept this quantity relatively constant. For this reason, 
they insisted that the valence CSV terms for up and down 
quarks in the nucleon be equal and opposite. Note that 
QCD evolution does not preserve this relation, but it is 
very nearly maintained for the region of evolution of in- 
terest. The overall coefficient, k, was then varied in a 
global fit to a wide range of high energy data. 

The value of k which minimised x 2 was k — —0.2. The 
MRST x 2 vs - K is shown as the bottom curve in Fig. [3] 
Clearly \ 2 nas a shallow minimum with the 90% confi- 
dence level obtained for the range —0.8 < k < +0.65. 
Since the form chosen by MRST for valence quark CSV 
is strongly constrained, and the resulting \ 2 minimum 
is quite shallow, one should not assign too much signifi- 
cance to their result. In global fits of this type, one should 
remember that, unless the shape of the constrained func- 
tion is in close agreement with the actual parton dis- 
tribution, the overall magnitude obtained in a global fit 
can be misleading. In Fig. 21 we plot the valence quark 
CSV PDFs corresponding to the MRST best fit value 
k = —0.2. They look extremely similar to the theoret- 
ical valence quark PDFs calculated by Rodionov et al. 
and shown in Fig. |2 they are also in good agreement 
with Sather's valence CSV distributions. This provides 
some theoretical support for the functional form chosen 
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Variation of j; with isospin violati 




FIG. 3: x profile for phenomenological isospin violating 
parton distributions, for sea quarks (top curve) and valence 
quarks (bottom curve), from the MRST group, Ref. Q. The 
quantity 5 associated with sea quark isospin violation is de- 
fined in Eq. 1101 while the coefficient k is defined in Eq. |§] 




FIG. 4: The valence quark CSV function from Ref. cor- 
responding to best fit value k, = —0.2 defined in Eq. Solid 
curve: xSd v (x); dashed curve: xSu v (x). 



by MRST. However, within the 90% confidence region 
for the global fit, the valence quark CSV PDFs could be 
either four times as large as that predicted by Sather or 
Rodionov, or it could be three times as big with the op- 
posite sign. The great value of the MRST global fit is 
that CSV distributions with this shape, and with values 
of k within this range, will not disagree seriously with 
any of the high energy data used to extract quark and 
gluon distribution functions. 

The MRST group also searched for the presence of 
charge symmetry violation in the sea quark sector. 
Again, they chose a specific form for sea quark CSV, 
dependent on a single parameter, i.e., 



u n {x) 
d n {x) 



<F{x) [1 + S] 
u p (x) [1 - 6} 



With the form chosen, the net momentum carried by 
antiquarks is approximately conserved; violation of mo- 
mentum conservation was found to be very small in the 
kinematic region of interest. 

Somewhat surprisingly, evidence for sea quark CSV 
in the global fit is substantially stronger than that for 
valence quark CSV. As shown in the top curve in Fig. 
13 the best fit is obtained for S = 0.08, corresponding 
to an 8% violation of charge symmetry in the nucleon 
sea. The x 2 corresponding to this value is substantially 
better than with no charge symmetry violation, primarily 
because of the improvement in the fit to the NMC fi — D 
DIS data HEl when u n is increased. The fit to the E605 
Drell-Yan data is also substantially improved by the 
sea quark CSV term. 

We can check the MRST results by comparing them 
with the charge ratio measurements that were summa- 
rized in Sect.[H] Eq. prelates the CSV parton distribu- 
tions to the charge ratio through 



Rc 



l 



3:r 



10Q(x) 



(5u(x) + 5u(x) — 5d(x) — 5d(x)j 



(10) 



We have taken the valence and sea quark PDFs from 
MRST j^], and in Fig. [5] we plot them against the exper- 
imental charge ratio obtained using the CCFR charged- 
current neutrino structure function and the NMC muon 
structure function. The long-dashed curve corresponds 
to the sea quark CSV term with the best value S = 0.08 
from the MRST fit. The remaining curves correspond 
to various values for valence quark CSV. The solid curve 
corresponds to the best value k = —0.2, the short-dashed 
curve corresponds to k — —0.8, and the dash-dot curve 
corresponds to k = +0.65. The latter two correspond to 
values of k at the 90% confidence level for the MRST fit 
to valence quark CSV. For values x < 0.085, we have not 
included the charge ratio experimental data since these 
changed significantly in the CCFR re-analysis. How- 
ever, the re-analyzed data resulted in a charge ratio that 
agreed with unity at about the 2-3% level, even down to 
x « 0.01. 

Both the valence and sea quark CSV distributions of 
MRST are in good agreement with the experimental lim- 
its from the CCFR/NMC data. At every x value the 
phenomenological values are within two standard devi- 
ations of the data. This is understandable, since both 
the NMC and CCFR data were included in the MRST 
global fit that extracted the CSV parameters. We have 
not included overall normalization errors, of roughly 2%, 
on both the CCFR and NMC data. Since the best fit 
valence quark PDFs of MRST (the solid curve in Fig. |3J) 
are very close to the theoretical predictions of Sather and 
Rodionov, the theoretical predictions are well within the 
upper limits on parton CSV from the best experiments. 
CSV effects in reasonable agreement with high energy 
data are substantially larger than predicted by theory; 
valence CSV effects could be four times as large as pre- 
dicted by Sather or Rodionov (or three times as large 
with the opposite sign), and sea quark CSV effects are 
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FIG. 5: The "charge ratio" R c — 1 defined in Eq. |21 com- 
pared with charge symmetry violating PDFs obtained by the 
MRST global fit, Ref. 0- Data is obtained from CCFR neu- 
trino cross sections of Ref. and NMC muon DIS,Ref. 0i. 
Long-dashed curve: sea quark CSV with best fit parameter 
S = 0.08 in Eq. I1UI solid curve: valence quark CSV corre- 
sponding to best fit value n = —0.2 in Eq. [5] short dashed 
curve: valence quark CSV for k — —0.8; dash-dot curve: va- 
lence quark CSV for k = +0.65. 



also significantly larger than theoretical predictions |40j . 

The only remaining issue is whether the theoretical 
and phenomenological PDFs agree with the limits on 
parton momentum, estimated by MRST as relative er- 
rors of about 4%. The best value obtained by MRST 
was K = —0.2. With this value, one obtains the relative 
momentum change for valence neutrons 



parton distribution functions, but in this ratio a large 
number of partonic corrections either cancel or are min- 
imized. 

The NuTeV group has measured neutrino charged- 
current and neutral-current cross sections on iron £| 
and extracted a value for sin 2 9 W equal to 0.2277 ± 
0.0013 (stat)±0.0009 (syst). This value is three standard 
deviations above the measured fit to other electroweak 
processes, sin 2 9 W = 0.2227 ±0.00037 42]. Such an effect 
can be interpreted as a 1.2% decrease in the left-handed 
coupling of light valence quarks to the weak neutral cur- 
rent. In addition to the Weinberg angle, several of the 
parameters of the Standard Model are constrained by the 
very precise measurements at electron-positron colliders 
|42| . Indeed, many of these parameters are determined at 
about the 0.1% level. Consequently, the discrepancy be- 
tween NuTeV and electromagnetic measurements of the 
Weinberg angle is surprisingly large. It suggests evidence 
of physics beyond the Standard Model, although David- 
son et ai, |43j have shown that it is quite difficult to pro- 
duce "new physics" (i.e., beyond the Standard Model) 
that fits the NuTeV experiment without violating other 
experimental constraints. In this paper, we will examine 
isospin- violating corrections to the NuTeV experiment. 

The NuTeV group did not directly measure the PW ra- 
tio, but instead measured the neutral-current to charged- 
current ratios R" and R v . These quantities, and their 
relation to the PW ratio, are given by 



2\SU V 



2(5 A, 



U y + D V U v + D v 



1% 



(11) 



r vN, 



This value is well within experimental limits. Even 
the largest value of K within the 90% confidence limit, 
K — —0.8, corresponds to a momentum change of 4%. 
Thus all of the phenomenological valence CSV distribu- 
tions (and the theoretical distributions of Sather and Ro- 
dionov et al.) correspond to momentum values within 
experimental limits. 



V. CHARGE SYMMETRY VIOLATION AND 
NEUTRINO DIS REACTIONS 



R~ 



(R u - i? 57 ) 
l-rR v 



R v = 



p 2 M n c°) 



(13) 



After measuring the ratios R v and R", the NuTeV group 
obtain the Weinberg angle through a detailed Monte 
Carlo simulation of the experiment. As a result, it is 
substantially more difficult to estimate the errors due to 
corrections to this exp eriment . For partonic corrections, 
the NuTeV group [44J have provided functionals that es- 
timate the change in an experimental quantity £ due to 
changes in a parton distribution 5(x): 



In 1973, Paschos and Wolfenstein |4lJ suggested that 
the ratio of neutral-current and charge-changing neutrino 
cross sections on isoscalar targets could provide an inde- 
pendent measurement of the Weinberg angle (sin 2 9 W ). 
The Paschos- Wolfenstein (PW) ratio R~ is given by 



R~ = 



P% 



(<Jcc°) 



1 • 2 

= sin 

2 



(12) 



In Eq.U2 (cr^c ) and i^cc ) are respectively the neutral- 
current and charged-current inclusive, total cross sections 
for neutrinos (or antineutrinos) on an isoscalar target. 
The quantity po = M W /(M Z cos9 w ) is one in the Stan- 
dard Model. Not only is the PW ratio independent of 



AS 



F [£, S; x] S(x) dx. 



(14) 



In Eq.^1 F [£, 5; x] is the functional, and the net change 
in the observable is obtained by integrating over Bjorken 
x. The NuTeV collaboration provided functionals appro- 
priate for both valence and sea quark CSV. 

The correction to the Paschos- Wolfenstein ratio arising 
from isospin violation in the parton distribution functions 
has the form 



AR C 



1 7 • 2, 

1 sin ( 

3 

SUy - 5D V 
2(U V + D v ) 



4a s / 1 . 2 

sin 9 W 

9tt \2 
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where 6Q V 



x Sq v (x) dx . 



(15) 



Note that the isospin violating correction to the PW ratio 
arises entirely from charge symmetry violation in the va- 
lence PDFs. Inserting t he q uark model results of Sather 
36] and Rodionov et al. jjtfl into En.HBI the CSV correc- 
tion to the PW ratio is found to be AR CSV = —0.0021 us- 
ing Sather's parton distributions, and AR CSV = —0.0020 
using the CSV distributions of Rodionov et al, |45| . 
These are theoretical predictions, as the Sather and Ro- 
dionov papers were published in 1992 and 1994, respec- 
tively, well before the analysis of the NuTeV data. Like 
the isoscalar corrections to the PW ratio, the CSV correc- 
tions involve the ratio of identical moments of parton dis- 
tributions, and are thus independent of Q 2 ; consequently 
they do not depend on any details of QCD evolution. 

The negative sign of the result means that CSV correc- 
tions will decrease the discrepancy between the NuTeV 
result for the Weinberg angle, and the Weinberg an- 
gle extracted from the extremely precise data obtained 
from electron-positron colliders near the Z mass. Indeed, 
the predicted CSV corrections to the PW ratio remove 
roughly 40% of the magnitude of the NuTeV anomaly. 
As discussed in the preceding Section, less than 1% of 
the valence quark momentum is carried by the relevant 
CSV combination, so this is well within the experimental 
limits on quark momentum. 

Londergan and Thomas [46( showed that Sather's for- 
mula of Eq. [7] provided an analytic estimate of CSV cor- 
rections to the PW ratio; taking the second moment of 
Sather's equations gives 

SD V = §D V+ 6 -^; SU v = §(U v -2) (16) 

Sather's approximation predicts that the CSV correction 
to the PW ratio depends only on the fraction of nucleon 
momentum carried by up and down valence quarks. Eq. 
EDp re dicts that SD V will be positive and SU V negative, in 
agreement with the theoretical model results. It is espe- 
cially interesting that Sather's analytic approximation al- 
lows one to calculate the CSV correction to the PW ratio 
directly from valence quark PDFs, without ever having 
to calculate specific CSV distributions. Using Eq.^3 wc 
also calculated the CSV correction to the PW ratio us- 
ing the CTEQ4LQ phenomenological parton distribution 
Hg at Q 2 = 0.49 GeV 2 , obtaining AR CSV = -0.0021 - 
again in agreement with the results obtained from the 
two quark model calculations. Eq. 1161 shows why the 
CSV corrections are so similar in various models; the 
correction depends only on the momentum carried by up 
and down valence quarks, a quantity well determined in 
PDFs. 

Since the NuTeV group 0, ^| did not directly mea- 
sure the PW ratio, one must multiply the CSV PDFs 
with the relevant functional in Eq. [21 44]. This re- 
quires evolving the CSV parton distribution functions 
from the quark-model scale to Q 2 — 20 GeV 2 , the aver- 



age momentum transfer for the NuTeV experiment. Af- 
ter QCD evolution and folding with the NuTeV function- 
als, the resulting CSV corrections to the NuTeV result 
are AR CSV = -0.0015, -.0017 and -0.0014 for the Ro- 
dionov, Sather and CTEQ4LQ PDFs, respectively. The 
CSV corrections decrease the NuTeV discrepancy in the 
Weinberg angle by about 30%, or one standard deviation. 
While some groups have obtained substantially smaller 
estimates for the charge symmetry violating contribution 
to the NuTeV anomaly |4J> |4g , in the first case the con- 
straint on baryon number was not respected while the 
second paper assumed a completely different mechanism 
for parton CSV. On the other hand, as we discuss next, 
the phenomenological values for CSV allow considerably 
larger isospin violation than any of the theoretical esti- 
mates. 

Valence quark CSV makes a substantially larger con- 
tribution than sea quark CSV to the extraction of the 
Weinberg angle from neutrino DIS. Using the sea-quark 
CSV and the best-fit value for valence quark CSV ob- 
tained by the MRST group, would remove roughly 1/3 
of the NuTeV anomaly. The value k = —0.6, within the 
90% confidence limit found by MRST, would completely 
remove the NuTeV anomaly, while the value k = +0.6 
would double the discrepancy. The MRST results show 
that isospin violating PDFs are able to completely re- 
move the NuTeV anomaly in the Weinberg angle, or to 
make it twice as large, without serious disagreement with 
any of the data used to extract quark and gluon PDFs. 

There are other possible corrections to the Weinberg 
angle extracted by the NuTeV experiment, but for the 
most part the corrections are small and/or well under- 
stood. The correction for the neutron excess in iron is 
large, but apparently well known [43l 0] . Radiative cor- 
rections are also large, but also believed to be under 
control although the radiative corrections are be- 
ing re-analyzed |50f. The situation with strange quarks 
is not certain. Spontaneously broken chiral SU(3) x 
SU(3 ) sy mmetry implies an asymmetry between s(x) and 
s(x) [5l| . As a direct consequence of such an asymmetry 
there is a correction to the extraction of the Weinberg 
angle that depends on the momentum asymmetry 



Sy 



x [s(x) — s(x)] dx . 



(17) 



A positive (negative) value for S v , meaning that strange 
quarks carry more (less) net momentum than strange 
antiquarks, would decrease (increase) the size of the 
Weinberg anomaly. The strange quark and antiquark 
distributions can be extracted from the cross sections 
for opposite sign dimuons in reactions induced by neu- 
trinos or antineutrinos; such reactions have been mea- 
sured by the CCFR |3 and NuTeV IX] collaborations. 
The NuTeV group has extracted strange and antistrange 
PDFs from the dimuon production reactions; they ob- 
tain values for S v that are consistent with zero, or per- 
haps small and ne gati ve [i^]. The CTEQ group on the 
other hand 153l " l54j . has included the dimuon data in 
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a global fit of PDFs; they obtain positive values for S v 
that could account for roughly one-third of the Weinberg 
angle anomaly. In both NuTeV and CTEQ analyses the 
strange quark distributions are determined from the same 
CCFR/NuTeV dimuon data, but at present the two anal- 
yses obtain qualitatively different results for the strange 
quark momentum asymmetry. 

In conclusion, it seems that the magnitude of CSV 
effects allowed by the MRST fit makes charge symme- 
try violation one of the only viable explanations for the 
anomalous value of the Weinberg angle obtained in the 
NuTeV neutrino experiment. If CSV effects are suffi- 
ciently large to remove the Weinberg angle anomaly, such 
effects should be visible in various other experiments. We 
had previously suggested several experiments that could 
potentially reveal the presence of parton CSV [jj. How- 
ever, the magnitude of those effects were based on the- 
oretical calculations that predicted substantially smaller 
CSV effects than are allowed by MRST. It is clearly of 



great interest now to investigate this issue experimen- 
tally. 
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